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Abstract. Hydrogen adsorption by graphite is examined by classical molecular 
dynamics simulation using a modified Brenner REBO potential. Such interactions 
are typical in chemical sputtering experiments, and knowledge of the fundamental 
behavior of hydrogen and graphene in coUisional conditions is essential for modeling 

Q ' the sputtering mechanism. The hydrogen adsorption rate is found to be dependent 

on the incident hydrogen energy and not on graphene temperature. Rather than 
destroying the graphene, hydrogen incidence at energies of less than 100 eV can be 

^ , classified into three regimes of adsorption, reflection and penetration through one 

CO ' or more graphene layers. Incidence at the lowest energies is shown to distort the 

\l , graphene structure. 
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j^ ■ 1. Introduction 

In chemical sputtering experiments, such as experiments using the diver tor of 
a nuclear fusion device, plasma-carbon interactions have been reported to yield 
hydrocarbon molecules of small number of atoms[l]. However, the mechanism of 
Q ' hydrocarbon generation has not been elucidated yet. In the present study, this 

O ' mechanism is examined in detail through simulation of the collision of hydrogen 

atoms with graphene as one of the fundamental processes in plasma-carbon inter- 
- action. 

rN \ A classical molecular dynamics (CMD) simulation scheme is adopted in the 

jrt ' present study to allow the dynamics of a many-particle system to be simulated 

over an adequate length of time given limited computer resources. Quantum me- 
chanical simulations could not be performed for the same time of interaction using 
the resources available. A new model potential between hydrogen and graphene is 
introduced for these simulations in order to incorporate Brenner's REBO potential 
[2] as the basis for the chemical reaction. The carbon wall is treated as graphene 
layers which are minimal structure of both a graphite and a carbon fiber. However, 
for simplification, they do not inter- layer interaction. The kinetic energy of incident 
hydrogen is set to less than 100 eV to facilitate comparison with the results of 
divertor experiments. 
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2. Simulation model 

The interaction potential in CMD has been developed through contributions by 
many theorists [2-5]. The Brenner's original REBO potential [2] has the following 
form: 

U=Y1 [y''in,)-h{{r})V^{n,)], (2.1) 

i,3{>i) 

where rij is the distance between atoms i and j, V is an attractive term, V^ 
is a repulsive term, and the function bij{{r^) includes all the effects of molecular 
orbitals. However, if chemical reactions occur, the REBO potential breaks energy 
conservation. To deal with chemical reactions, new functions expressing conjugation 
effects are thus proposed: 

carbon carbon 

Ntr' = l+ Y. nnu)CN{Nl)+ Y. nr,i)Cr,{Nl), (2.2) 

where f^ is a cut-off function for the distance between atoms, and 

( 1 ffx sC 2, 

Cn{x)^< [l + cos(7r(a;-2))]/2 ff2<a;^3, (2.3) 

[ ffa;>3. 



conj 



In contrast to Brenner's original formulation [2] , the second and third terms of N, 
are not squared. The tricubic spline functions F and T in Ref. 2 are redefined, as the 
original functions have N^°^-' as a variable. Denoting a new function of i^ by F', we 
obtain F'(z,j,4) = F(i,j,6), F'{i,j,k ^ 5) = F(i,j,9) and F'{i,j,k) = F{i,j,k) 
for the other cases. The spline T{i,j,k) in Ref. 2 is redefined in a similar way. 
The above modifications of N^""'-' , F' and T' provide differentiability at the cut-off 
point, conserving the total energy in the chemical reaction. This modified potential 
is employed in the present CMD simulations. 

Graphenes are set parallel to the x — y plane in z < region. Each graphene is 
comprised of 160 carbon atoms with a periodic boundary condition. The velocities 
of all atoms in the graphene are distributed according to a Maxwell's distribution at 
the graphene temperature. Incident hydrogen atoms arrive normal to the graphenes 
from the z — 4 A. The simulation continues until the incident atom is reflected, 
the atom leaves the graphenes into the z > region, or until the atom becomes 
trapped by the graphenes. 

A total of 200 simulations was performed for each graphene temperature and 
incident energy, with the initial position of the hydrogen atom on the z = 4 A plane 
varied randomly to obtain statistical information. The present CMD simulation 
satisfies the constant- volume-energy (NVE) condition. 



3. Results and Discussion 

In the present simulations of incident hydrogen energies Ei of less than 100 eV, the 
hydrogen atoms were found to be absorbed (or reflected) by the graphenes without 
destroying the graphenes. As shown in Fig. 1(a), the simulations suggest that the 
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Figure 1. (a) Hydrogen adsorption rate vs. graphene temperature, and (b) incidence 

energy at 800 K 

(a) 




o 


Hydrogen 


• 


Carbon 


— 


Bond 


= 


C overhang (u) 


^ 


CH distance (w) 



■1.0 -0.5 0.0 0.5 1.0 

u : C overhang [X] 



Figure 2. (a) Schematic of a hydrogen atom fixed above a carbon atom showing the over- 
hang from flat graphene (x) and the distance between the carbon atom and the incident 
hydrogen atom (j/). (b) Schematic of reaction paths in the hydrogen-graphene colUsion 
process in total potential energy space. The white regions denote higher potential en- 
ergy than eV. Adsorption, reflection and penetration are denoted by arrows 1, 2 and 3, 
respectively. 

hydrogen adsorption rate on graphene is independent of the graphene temperature. 
This result is not consistent with experimental facts[6, 7]. However, the tempera- 
ture defined in this simulation (T^) is not the same as the temperature measured in 
experiments (Te) due to the omission of inter-layer interaction in the present sim- 
ulation. The graphenes form a domain structure by the inter-layer interaction in 
the carbon wall. The experimental temperature T^ therefore includes all the kinetic 
energy of the relative motion of each domain, each layer, and each atom, whereas 
Ts accounts solely for the kinetic energy of atoms. 

Figure 1(b) shows that the hydrogen adsorption rate is dependent on Ei. The 
mechanism of adsorption can be classified into the following three types of processes: 

Adsorption at E^i < 5 eV: As the incident hydrogen atom approaches the 
graphene, the nearest carbon atom overhangs spontaneously from the flat graphene 
(Fig. 2(a)). Because, this deformation of the graphene leads the total potential 
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energy to the minimum point (arrow 1, Fig. 2(b)). The kinetic energy converted 
from potential energy is diffused into the entire atoms in the graphene through 
the overhanging carbon atom. The incident atom thus loses kinetic energy and is 
trapped. 

Reflection at 5 eV < _Ei < 50 eV: The velocity of the incident atom is so high 
under high-energy conditions that spontaneous overhanging of the carbon atom 
can not occur, and the total potential does not pass through the minimum point. 
The incident atom bounces off this potential barrier (arrow 2, Fig. 2(b)) and the 
incident atom is reflected. The adsorption rate under such conditions is therefore 
low. 

Penetration at Ei > 50 eV: When the energy Ei is greater than the poten- 
tial barrier, the incident atom penetrates the graphene (arrow 3, Fig. 2(b)). The 
incident atom gradually loses kinetic energy by repulsion of the carbon atom as it 
penetrating several graphene layers. The reaction thus eventually changes into a 
reflection process. After repeated bouncing between graphene layers, the incidence 
atom is trapped. 



4. Conclusions 

The collision processes between incident hydrogen and graphene at energies of less 
than 100 eV were investigated by CMD simulations using a modified REBO poten- 
tial. It was found that hydrogen does not destroy the graphene at these energies, but 
is instead absorbed or reflected by the graphene. The rate of hydrogen adsorption 
was shown to be dependent on the incident energy and not on the graphene tem- 
perature. The collision behavior was classified into the three processes; absorption 
of hydrogen in graphene (Ei < 5 eV), reflection by graphene (5 eV < £"1 < 50 eV), 
and penetration through one or more graphenes (£"1 > 50 eV). 
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